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The modified Chaplygin gas (MCG) is a strong candidate for the unified model of dark matter 
and dark energy. The equation of state of this modified model is valid from the radiation era to 
the ACDM model. In early epoch (when p was large), dark matter had the dominant role while 
at later stages (when p is small), the MCG model behaves as dark energy. In this work, we have 
found exact solution of static spherically symmetric Einstein equations describing a wormhole for an 
inhomogeneous distribution of modified Chaplygin gas. For existence of wormhole solution, there 
are some restrictions relating the parameters in the equation of state for MCG and the throat radius 
of the wormhole. Physical properties and characteristics of these modified Chaplygin wormholes are 
analyzed in details. 
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I. INTRODUCTION 

Recent observational evidences of anisotropy in the 
Cosmic Microwave Background radiation [1] and the 
data from type la Supernovae SN 1997H redshift [2-4] 
suggest that the universe is flat and is undergoing at 
present an accelerating phase preceeded by a period 
of deceleration. These observational facts can not be 
explained by the known ordinary (baryonic) matter 
or radiation. There should be a significant amount 
of energy density of the universe which should be an 
extraordinary non-baryonic matter (dark matter) and 
energy (dark energy). Also, observational facts indi- 
cate that, dark matter should be of the order of 25% 
of the critical density and dark energy contributes 
about I of the the critical density [4-7]. Further, 
the universe in very recent past (at redshift z ^ 1), 
is dominated by uniformly distributed dark energy 
with negative pressure and is probably responsible 
for the present phase of acceleration of the universe. 
In last few years, several models have been suggested 
to incorporate the recent observational evidences and 
among them, a single component perfect fluid having 
exotic equation of state, known as Chaplygin gas 
[8-18] is of great interest. A further generalized form 
is known as modified Chaplygin gas [the equation of 
state is given in Section III, equation (12)]. The pa- 
rameters A, B and a are universal positive constants. 
The special choice A — and a = 1 corresponds 
to Chaplygin gas, first introduced to describe lifting 
forces on a plane wing in aerodynamics process [12]. 
Its generalization with the choice A = 0, a > is 
known as generalized Chaplygin gas (GCG), first 
introduced by Kamenshchik et al.[13] and Bento et 
al.[14]. For small energy density, the GCG model will 
have negative pressure while at high energy density, 
the model behaves as an almost pressureless fluid. 
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However, MCG model, an extension of GCG model 
[8, 15], can interpolate states between standard fluids 
at high density and constant negative pressure at low 
energy. In fact, the MCG model for 7 = ^ describes 
the evolution from radiation epoch at early time to 
ACDM era at late time (where the fluid has constant 
energy density and behaves as a cosmological con- 
stant). This causes the present accelerating phase 
of the universe [8, 9, 15]. Another point that also 
goes against the GCG model is the measurement of 
the velocity of sound. In GCG model, the velocity of 
sound is negligible at early times and approaches the 
speed of light at late times. But in cosmic evolution, 
sound velocity comparable to the light velocity is not 
compatible. 

From phcnomcnological view point, MCG model is 
interesting and can be motivated by the brane world 
interpretation [14]. Further, this model is consistent 
with various classes of cosmological tests viz. gravi- 
tational lensing [17], gamma-ray bursts [18] as well 
as the above mentioned observations. Moreover, the 
present MCG model (also GCG model) is naturally 
constrained through cosmological observables [16]. 
Also, for low energy density, the present model has 
similarity with GCG model and so the equation of 
state is that of a polytropic gas [15] with negative 
index. Therefore, it is possible to have astrophysical 
implications of the present model with an alternative 
way of restricting the parameters [16]. 

As in MCG model, though the pressure is initially 
positive, but with the evolution of the universe, the 
pressure becomes negative and finally behaves as 
phantom energy. It is well known that phantom 
energy violates null energy condition, a natural 
scenario for the existence of traversable wormholes 
[19, 20]. This purely theoretical ingredient pos- 
sesses a peculiar property viz. "exotic matter", 
whose stress-energy tensor violates the null energy 
condition. This common fact of violation of null 
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energy condition motivates to study wormholes 

supported by MCG equation of state. In this work, 
we shall construct static and spherically symmetric 
traversable wormhole geometries, satisfying the MCG 
equation of state. We shall consider matching of 
these wormhole geometries to an exterior vacuum 
space time and examine physical properties and 
characteristics of these solutions. We shall also study 
the traversability conditions [19, 21] and measure the 
amount of averaged null energy condition violating 
matter [22, 23] in particular cases. 



vacuum solution. 

The following properties need to be imposed for the 
existence of a wormhole solution: [19, 20] 



I. As ro is the throat radius, so b{ro) = ro (2) 



II. A flaring out condition of the throat, i.e 



The remainder of our paper is organized as follows: 
Section II provides basic equations for wormhole 
geometry. In Section III, we have formulated the 
fundamental equations describing wormhole in 
modified Chaplygin gas and the restrictions on the 
parameters of the equation of state of MCG have 
been presented. Section IV deals with asymptotic 
flatness of the wormhole geometry. The non-flat 
(aymptotically) wormhole solutions are bounded by 
an exterior vacuum space time for which, details of 
surface stresses are shown. In section V, some specific 
wormhole solutions are derived and their physical 
properties and characteristics are explored. We have 
considered the choices viz. (i) a particular relation 
between $ and b (in differential form), (ii) constant 
redshift function, (iii) two typical choices of the 
shape function and lastly (iv) isotropic pressure with 
a typical form of 6(r). The traversability condition 
has been studied in details for constant redshift 
function with a toy example. Also for the linear 
form of 6(r), the amount of averaged null energy 
condition violating matter has been estimated using 
"volume integral quantifier". Finally, the paper ends 
with discussion and concluding remarks in Section VI. 



(6 - b'r) 



> 0, which at the throat simplifies to b'{ro) < 1 

(3) 



III. 



6(r) < r for r > ro 
(regularity of the metric coeflacient) 



(4) 



IV. Traversability criteria: For traversable worm- 
hole, there should not be any horizons present i.e, $ 
must be finite everywhere. 

Using the Einstein field equation G^u = SttT^,^, in 
an orthonormal reference frame (G = C = 1), we 
obtain the following set of equations 



p{r) 



b' 



87rr2 



Pr{r) = - 



r V r I r-' 



(5) 



(6) 



II. BASIC EQUATIONS FOR 
WORMHOLE GEOMETRY 

A spherically symmetric and static wormhole is 
given in Schwarzschild coordinates by the space time 
metric 



and 



1 

2r2 ' 



(7) 



ds^ = -exp [2$(r)] 



dr^ 



1 



b{r) 



(1) 



where the two arbitrary functions of r viz. $(r) 

and b{r) are as usual termed as the redshift function 
and the shape function respectively. The radial 
coordinate r ranges over [ro,oo) where the minimum 
value To corresponds to the radius of the throat of 
the wormhole. But there may be a cut-off of the 
stress-energy tensor at any finite junction radius R, 
where the interior space time matches to an exterior 



in which p(r) , Pr (r) and pt (r) are the energy density, 
radial pressure and tangential pressure respectively. 
The energy conservation equation gives 



= -{pt-Pr) - {p + PrW 



(8) 



This equation can be interpreted as the hydrostatic 
equilibrium equation for the material threading the 
wormhole. 
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The above Einstein equations can be rearranged to 
obtain metric coefficients in terms of the components 
of the stress-energy tensors as 



b + Snprr^ 
2r^ (l - ^) 



(9) 



(10) 



The above flaring out condition impHes that the 
wormhole should be threaded with matter violating 
the null energy condition (NEC) (a fundamental prop- 
erty of wormholes). Matter that violates the NEC is 
termed as exotic matter [19, 24, 25]. The mathemati- 
cal form of the null energy condition is Tf^,yK'^K'^ > 
for any null vector K^^. However, in an orthonormal 
reference frame, we may choose Rf^ = (1,1,0,0) so 
that the violation of NEC simplifies to p + pr- < i.e, 



1 

Stt 



b'r~b 



r I r 



< 



(11) 



Note that, if we evaluate at the throat considering 
the flaring out condition and the finite character 
of $(r), the above inequality (11) is automatically 
satisfied, i.e, the null energy condition is violated. 



III. MODIFIED CHAPLYGIN GAS: THE 
EQUATION OF STATE 

The MCG is described by the equation of state [8, 9] 



Pr 



< a < 1 



(12) 



where the parameters A and B are two universal 
positive constants. This equation of state shows a 
radiation era (when A = 1/3) at one extreme (when 
the matter density p is very large) and a ACDM 
model at the other extreme (when p is infinitesimally 
small). However, at all stages it shows a mixture. 
Further, somewhere in between the two extreme, 
there is one stage when the pressure vanishes and 
the matter content is equivalent to pure dust (this 
cosmological model can be considered from the field 
theoretical point of view by introducing a scalar field 
having self- interacting potential [9]). 

Now eliminating p and Pr among the equations (9), 
(10) and (12), wc have the relation between $ and b 
(in differential form) as 



b + Arb' 



(b'Y 



2r2 (1 



(13) 



Any solution of the metric (1) satisfying relation 
(13) is termed as modified Chaplygin wormhole 
provided it obeys the restrictions (mentioned in the 
previous section) for wormholes. In principle, we 
now have four equations (8)-(10), (13) containing 
five unknowns viz. two metric coefficients ^(r) and 
b{r) and three physical quantities (components of 
stress-energy tensor) p, Pr and pt- So we may restrict 
one of the above unknown quantity (or a relation 
among them) for a specific solution. In the present 
problem, wc shall choose b{r) or ^(r) so that the 
solution fulfills the criteria of wormholes. Note that, 
as we shall consider the matter density p to be 
positive, so from equation (5), the shape function 
b{r) should be such that b'{r) > V r. 

Now at the throat r = ro, the components of the 
stress-energy tensor are (for the choice a = 1) 



(sTrro^) 



r(o) 



_ l-b'(ro) J_\ 



(14) 

The restriction 5'(ro) < 1 for the wormhole geome- 
try restricts the parameter A and B in MCG equation 
of state as 



B < 



(87rrg)2 



(15) 



In the present context, the violation of null energy 
condition i.e, p + Pr < demands the energy density 
to be restricted as 



B 



1 + A 



< 



(16) 



This is not unusual in MCG cosmology models. 
Usually in FRW cosmology, the evolution of p is given 
by (see eqn. (8) in Ref. [9]) 



B 



C 



l + A a3(i+^)(i+«) 



(l/l+a) 



(17) 



where a is the scale factor in FRW model and C is 
arbitrary integration constant. So if C is chosen to 
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be positive, then p"^^ > i.e, dominant energy 

condition is satisfied. But for the choice of negative 
values of C, the restriction (16) is automatically sat- 
isfied, a fundamental ingredient in wormhole physics. 
Moreover, the expression for the velocity of sound is 



dp 



A + 



Ba 



(18) 



The condition that Vg should not exceed the velocity 
of light demands 



then one must have 

R > Tb = 2M (event horizon). 

Using the junction conditions due to Darmois-Israel 
formalism [27], the components of the surface stresses 
of a dynamic thin shell [28, 29] are 



a (surface energy density) = 



4ttR 



2M 



Ba 

1-A 



(l/l+a) 



(19) 



Thus combining (16) and (19) we have the con- 
straint 



Ba 
1-A 



< p"+^ < 



B 



1 + A 



< 



(20) 



It should be mentioned that in the domain of 
exotic matter, |2 may not be interpreted as the 
speed of sound due to the lack of knowledge of a 
microphysical model of exotic matter. On the other 
hand, there are examples of exotic behaviour such as 
Casimir effect and the false vacuum [26] where < 0. 



IV. ASYMPTOTIC FLATNESS AND 
STRESS-ENERGY TENSOR CUT-OFF 

In the above, the basic structure of the solution has 
been shown and the restrictions for the parameters 
involved are presented. A solution for which $ ^ 
and ^ ^ as r ^ 00, corresponds to asymptotically 
flat space times. However, every solution may not 
be related to asymptotically flat space time. Then 
the stress-energy may be cut-off at a finite radial 
coordinate R and the interior solution of metric (1) 
is matched to an exterior vacuum space time at the 
junction interface r = R. The junction r = R may 
be a thin shell (if there are surface stresses at the 
junction surface) or simply a boundary surface (zero 
surface stresses). Hence for solutions corresponding 
to asymptotically non-flat geometry, the MCG 
distribution is restricted to the neighbourhood of the 
throat of the corresponding wormhole solution. 



R 



(22) 



and 



V (tangential surface pressure) 



87ri? 



1 - 2M ^ ^2 



(1 + R^'){1 -^-i^ + R^)+RR- ^^j^-y 



(23) 

The above expressions for the components of 
stress-energy tensor can be simplified to a great 
extent for the static case [30] by imposing d = = a, 
where the overdot denotes derivative with respect to 
the proper time r. 

Further, if Mg = 47rr^cr denotes the surface mass 
of the thin shell (for static case), then the total mass 
Mt of the wormhole has the expression [21] 



biR) 



2R 



(24) 



Also for the static case, if we have the boundary 
surface i.e, a = = V, then one obtains the following 
relationships 



b{R) = 2M and R = 2M [ ^-r^ 1 (25) 



As i? > 2M, so the redshift parameter C 
[=1+ R^'{R)] is restricted by C > 1- 



If the exterior vacuum space time is chosen for sim- 
plicity as the Schwarzschild space time i.e. 



V. SPECIFIC WORMHOLE SOLUTION 



dsl = -{l- 



2M\ 



1 , 9 df"^ 



-r^rffia^ (21) 



A. A special relation between the redshift 
function and the shape function: 
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fig. (a) shows the variation of 6(r) (given by eq.(27)) over 
the radial coordinate r for /it = 1 and ro = | . 



Suppose we choose 



violated. However, for traversability criteria for con- 
stant $, wc need only b'{r) at the throat, which can 
be obtained from (29) as (for a = 1) 



b'iro) = 



a/I + 12 Afxr* - 1 



2A 



(30) 



Note that for the flaring out condition we must have 



<1 + A, n = 



(31) 



We shall now examine an important concept of 
wormhole physics namely the traversability conditions 
for a human being to go forward through the worm- 
hole. The basic criteria for a smooth journey is that 
the acceleration felt by the traveler should not exceed 
earth's gravity ge (for details see ref. [19-21]) i.e. 



A b' 



i.e, $ 



/ 



A b' 



dr (26) 



2 (r - 6) ' J 2 (r-b) 
Then from equation (13), the solution for b becomes 



with 3/xr^ < 1 for 6'(ro) < 1. 



(27) 



Clearly, this solution does not correspond to an 
asymptotically flat space time. So the MCG is con- 
fined in a neighbourhood of the throat and the above 
solution can be matched to exterior vacuum geometry 
(as described in Section IV). Further, the restriction 
b(r) < r yields a range of r (for which wormhole solu- 
tion is possible) viz. 



ro < r < 



4/i 



1/2 



(28) 



and is shown diagramatically in fig. (a). 

B. Constcint redshift function: *'(r) = 

This choice of the redshift function gives (from (13)) 
the differential equation in the shape function as 



One can immediately inspect that b{r) 



is a particular solution for b{r) but it is not a worm- 
hole solution as b'{ro) < 1 (fiaring out condition) is 



*(7e*)' 



< Qe 



(32) 



Moreover, for smooth journey across the worm- 
hole, the tidal acceleration also should not exceed the 
earth's gravitational acceleration (gg) and we have the 
following restrictions [19] 



1-'- 
r 



7 

2r2 



b' 



^,2 _ ib^r-bl^, 



2r{r - b) 



+ 2{r - 6)$' 



\v' I < 9e (33) 



I < 9e (34) 



Here the separation between two arbitrary parts of 

-/ 

the traveler's body is denoted by jry* | and is chosen 
to be 2m along any spatial direction in the traveler's 
reference frame as in the literature [21, 25]. 

Suppose a human being travels with velocity v along 
the radial direction such that u > in —li < I < I2 
but w = at Z = — Zi and I2, where the two space 
stations arc assumed to be located (just outside the 
junction shell radius R) at proper radial distances I = 
—li and I = I2 respectively. {I = J ^ is the proper 



radial distance). Then the time of travel across the 
wormhole is given by [19, 21, 31] 



b{bT + Ar{b'r+^ - {8'K)"+'^Br^"+^ = (29) (AT), 



-h 



and (Ar)o = 



dl 



(35) 



where 7 = [1 — v'^]~^^'^ and (AT)t is the proper 
time interval of the journey recorded by the traveler's 
clock, (AT)o is the coordinate time interval recorded 
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by observers situated at the stations. 

For simplicity, if we choose wormholes correspond- 
ing to $ =constant (i.e, constant redshift), then for 
a traveler moving with non-relativistic velocity (i.e, 
7 « 1), the restrictions (32) and (33) are identically 
satisfied while the inequality (34), when evaluated at 
the throat, gives the upper limit of the velocity as 



■29e 



[1 - b'iro)] 



(36) 



As a toy example if we choose b'{ro) = ^, ro = 
lO^m, then one obtains w w 4 x 10^ m/s (choosing 
the equality sign) for the maximum traversal velocity. 
Further, choosing the junction radius R ~ lO^'m, 
the traversable time is (AT)t ~ (AT)o ~ ^ w 50 sec. 

C. Specific choice for tiie siiape function: 



I. b{r) = roi^y 



The particular choice of the shape function 



as there is an event horizon at the throat r = ro- 

However, it is possible to make the wormhole solution 
traversable by imposing the restriction 



1 + An) _ 32Tr'^B4 
2(1 - n) ~ n(l - n) 

Then $(r) simplifies to 



(40) 



$(r) = $0 - 



1 + An) 
2(1 - n) 





r 


In 






ro 



az 



a(a - 1) z"^ 



(41) 



As <l>(r) does not approach to a constant limit 
as r oo, so the solution is not asymptotically 
flat. Hence we may match this interior solution to 
an exterior vacuum space time at a junction radius 
R and $o can be determined from the boundary 
conditions of $ at the junction interface. 



b{r) = ro , < n< 1 

satisfies the wormhole criteria viz. 



(i) 6(ro) = ro , {ii) b'{r) =n[ — ] i.e, 



(37) II- b(r) =ro + d{r - rp) 



This typical choice of b{r) satisfies as above all the 
criteria to describe a wormhole provided < d < 1. 
It is clear that - — > as r ^ oo and the solution for 

r 

$ becomes 



b'{ro)=n<l, (m) 1-^ = 1- f^)' "<1 

r \ r J 

(38) 



Also^ = (^)^ 



as r — > OO. 



Then solving equation (13) we have 



^, , ^ [l + Ad)^, , 1 327r2B 



4 ln|r -ro\ + - (r - Tq)"* + ■^rl{r - ro) 



$(r) = $0 



{l + An) 
2(1 -n) 



In 



In 



l-n 



'S2n^Br^ 
n(l — n) 



In 



- 1 -1 

ro 



+ az 



+ 



a{a — 1) z"^ 



(39) 



where ^ = [^) — 1, a = and $o is the 
integration constant. It can be seen easily that the 
above solution describes a non traversable wormhole 



+3ro(r-ro) + - ro{r - ro) 



(42) 



which clearly describes a non-traversable wormhole. 
However, the choice 

A = i [64n^Br^ - l] 

makes $ finite at r = ro and the solution describes 
a traversable wormhole (but not asymptotically flat). 
So the MCG is conflned around the throat up to a 
junction interface r = R. 

We shall now examine how much of matter is in the 
space time, which violates the averaged null energy 



condition. According to Visser et al.[22, 23], this in- 
formation can be obtained by "volume integral quan- 
tifier" defined as 



\n\{A+l)pl-B\-\\n 



^ — I 2(1+A) "'IV-- ' -;ru ^\ 2 

Suppose we choose the shape function as [32] 



_(l±i)£L 
A+1)pI-B 



W)+Pr{r)] dV 



I 

J rn 



(r-6) 



In 



^2$ 



1 - 



dr 



(43) 



Now using the solution (42) for $ and the choice 
6(r) = ro -|- d{r — ro), the above integral simplifies to 



= -(l + Ad) 



,2 a 



+ [l + Arf-r^(l-rf)]a 



(44) 



where a = R — ro is the radial distance of 
the junction interface from the throat. Thus as 
a ^ 0, ly ^ 0. This shows that, similar to the 
phantom [21] or GCG wormholes [25], it is possible 
to construct a wormhole solution in MCG, where the 
amount of modified Chaplygin gas is arbitrarily small. 
This conclusion is in agreement with the results in 
the literature [21, 25] regarding Chaplygin gas. The 
interesting point to be noted is that, theoretically 
it is possible to construct wormholes supported by 
infinitesimal amount of exotic fluids in cosmology, 
which may responsible for the present accelerating 
phase of the universe. 

D. Isotropic Pressure: 



In case of isotropic pressure (i.e, pr = pt), the en- 
ergy conservation equation (8) simplifies to 



(Ap"+i + Ba) 



[(A+l)/j"+i 
which on integration gives 



B]p 



(45) 



$ = - 



A 



(A + l)(a + l) 



\n\{A + l)p 



a + l 



B 



+ 



a 



a + l 



In 



{A + l)p"+i 



{A + l)p«+i - B 



+ c 



(46) 



The constant c can be determined from the fact 
that, at the throat r = ro, 



PO = 2A 



(^)2+AAB-^ [see eqn.(14)] 



and (j){r = ro) = (po (say), i.e, (for a = 1) 



b{r) = ro + 7^ro 1 - — 
L r . 

with < 7^ < 1. 

Then from field equation (5) we have 



(47) 



Pi"-) = ^ 



pi = 



Stt 



(48) 



and the explicit expression for $ becomes 



$(r) = c - 



(^ + l)(a + l) 



In 



a+l 



j.1(pt+l) 



B 



a+ 1 



In 



a+l 



(A+ l)/j!;'+i -r4("+i)B 



(49) 



Clearly, the solution is not asymptotically flat and 
the radius of the junction interface should be re- 
stricted by 



R< 



B 



(50) 



Finally, it should be mentioned that wormhole 
solution for constant energy density is identical to 
that in ref. [25], so we have not presented it here. 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

The modified Chaplygin gas model is a unified 
model of dark matter and dark energy, which 
causes the present accelerated expansion of the 
universe. First of all in the context of cosmology, 
the MCG model describes a smooth transition from 
a decelerated expansion to the present epoch of 
cosmic acceleration (i.e, a continuous transition 
from radiation era to ACDM model). Also from the 
phenomenological view point, the model describes 
unified macroscopic prescription of dark matter and 
dark energy. Lastly, the model is favorable than 
the GCG model (clearly presented in the introduc- 
tion) and is consistent with several cosmological tests. 

In the present work, we have examined whether 
wormhole solutions (traversable or non-traversable) 
are possible in the background of MCG model. It is 
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found that, for the existence of a generic solution of 
a wormhole, the parameters in the equation of state 
of MCG should be restricted as (for a = 1) 

and consequently, there is violation of null energy 
condition (a fundamental requirement in wormhole 
physics). We have presented four specific worm- 
hole solutions and analyzed the physical properties 
with characterization of these modified Chaplygin 
wormholes. In the first solution, we have assumed a 
differential relation between the redshift function <l>(r) 
and the shape function b{r). The wormhole solution 
is valid within a finite range of r which is matched 
with an exterior vacuum space time. The range of r 
for the wormhole solution is shown graphically. For 
the choice of constant redshift function, no analytic 
form of the shape function is possible. Only we have 
determined the restriction in the parameters involved 
to satisfy the traversability criteria. Then solutions 
are obtained for two specific choices of the shape 
function. In both cases, the wormhole solutions can 
be made traversable by imposing certain restriction 
on the parameters involved. For the second choice, 
we have calculated 'volume integral quantifier' (/„) 
to have an idea about the amount of matter violating 
averaged null energy condition. We have seen that 
^ as the junction interface approaches to the 



throat. Lastly, we have obtained wormhole solution 

for isotropic pressure. Here also due to complicated 
form of $ a function of energy density p, the form 
of the shape function is assumed and solution is 
confined in a finite restricted region. Finally, we 
have shown in general the traversability condition 
for a human being. The restrictions are simplified 
a great deal for constant redshift. We have finally 
estimated the time of travel for a traveler moving 
with non-relativistic speed. It should be mentioned 
that the space stations should be located at large r, 
so that the geometry at the stations is nearly flat [31]. 
As a concluding remark, we may say that the MCG 
traversable wormholes have far-reaching physical 
and cosmological implications. In addition, for the 
use of interstellar shortcuts, an absurdly advanced 
civilization may use them as time machines (probably 
violating causality) [20, 33, 34] (for review see [35]). 
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